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CONSERVATION LAWS OF RANDOM MATRIX THEORY 

NICHOLAS M. ERCOLANI * 

Abstract . This paper presents an overview of the derivation and significance of recently derived 
conservation laws for the matrix moments of Hermitean random matrices with dominant exponential 
weights that may be either even or odd. This is based on a detailed asymptotic analysis of the 
partition function for these unitary ensembles and their scaling limits. As a particular application 
we derive closed form expressions for the coefficients of the genus expansion for the associated free 
energy in a particular class of dominant even weights. These coefficients are generating functions for 
enumerating g-maps, related to graphical combinatorics on Riemann surfaces. This generalizes and 
resolves a 30+ year old conjecture in the physics literature related to quantum gravity. 

o 

1. Introduction 

This paper will present an overview of some recent developments in the appli- 
cation of random matrix analysis to the topological combinatorics of surfaces. Such 
applications have a long history about which we should say a few words at the outset. 
The combinatorial objects of interest here are maps. A map is an embedding of a 
graph into a compact, oriented and connected surface X with the requirement that 
the complement of the graph in X should be a disjoint union of simply connected 
^^ open sets. If the genus of X is g, this object is referred to as a g-map. The notion of 

g-maps was introduced by Tutte and his collaborators in the '60s [37] as part of their 
C* investigations of the four color conjecture. 

In the early '80s Bessis, Itzykson and Zuber, a group of physicists studying 't 
Hooft's diagrammatic approaches to large N expansions in quantum field theory, dis- 
covered a profound connection between the problem of enumerating g-maps and ran- 
dom matrix theory [4|. That seminal work was the basis for bringing asymptotic 
analytical methods into the study of maps and other related combinatorial problems. 
Subsequently, in the early '90s, other physicists [HI E2] realized that the matrix 
_i model diagrammatics described in [4] provide a natural means for discretizing the 

Einstein-Hilbert action in two dimensions. From that and a formal double scaling 
VsQ limit, they were able to put forward a candidate for so-called 2D Quantum Gravity. 

~f-T This generated a great deal of interest in the emerging field of string theory. We refer 

to [10] for a systematic review of this activity and to [30] for a description of more 
C\| recent developments related to topological string theory. 

All of these applications were based on the postulated existence of a 1/n 2 asymp- 
^- totic expansion of the free energy associated to the random matrix partition function, 

where n denotes the size of the matrix, as n becomes large. The combinatorial sig- 
nificance of this expansion is that the coefficient of 1/n 29 should be the generating 
function for the enumeration of g-maps (ordered by the cardinality of the map's 
vertices) . In |14| the existence of this asymptotic expansion and several of its impor- 
tant analytical properties were rigorously established. This analysis was based on a 
Ricmann-Hilbcrt problem originally introduced by Fokas, Its and Kitaev [20] to study 
the 2D gravity problem. 

The aim of this paper is to outline how the results of [H] and its sequel [16j have 
been used to gain new insights into the map enumeration problem. In particular, 
we will be able to prove and significantly extend a conjecture made in [1] about the 
closed form structure of the generating functions for map enumeration. 



X 



'Department of Mathematics, The University of Arizona, Tucson, AZ 85721-0089, 
(ercolaniSmath.arizona.edu). Supported by NSF grant DMS-0808059. 
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Over time combinatorialists have made novel use of many tools from analysis 
including contour integrals and differential equations. In this work we also introduce 
nonlinear partial diferential equations, in particular a hierarchy of conservation laws 



reminiscent of the shallow water wave equations [3 8) (see (6.7|). This appears to make 



contact with the class of differential posets introduced by Stanley [35j (see Remark 



4.9) 



2. Background The general class of matrix ensembles we analyze has proba- 
bility measures of the form 



1 



d»ti = ^ul I ^ M Tr[^(M,i,)] \dM, where (2.1) 



V j (X;t j ) = h 2 + t f^ (2.2) 

defined on the space T-L n of n x n Hermitcan matrices, M, and with g s a positive 
parameter, referred to as the string coefficient. The normalization factor Z^ n '{g s ,tj), 
which serves to make fit & probability measure, is called the partition function of this 
unitary ensemble. 

Remark 2.1. In previous treatments, \1^\ \16\ WA \17j , we have used the parameter 
\/N instead of g s . This was in keeping with notational usages in some areas of 
random matrix theory; however, since here we are trying to make a connection to 
some applications in quotum gravity, we have adopted the notation traditionally used 
in that context. This also is why we have scaled the time parameter tj by 1/j in this 
paper. 

For general polynomial weights V it is possible to establish the following funda- 
mental asymptotic expansion [T3] , [TS] of the logarithm of the free energy associated 
to the partition function. More precisely, those papers consider weights of the form 

VW^ + Jl-^- (2-3) 

e—i 

with j even. 

We introduce a renormalized partition function, which we refer to as a tau function 

representation, 

r 2 (t)= Z( " )(f) (2 4) 

where t= (t±,...tj) £W . The principal object of interest is the large n asymptotic 
expansion of this representation for which one has the result [141 I16j 



logr 2 (f) = n 2 e (x,t) + e 1 (x,t) + -^e 2 (x,t) + ■■■ + -n—~e g {x,t) + ... (2.5) 
' J n n 9 

as n— ^co while g s — >0 with x = ng s , called the 't Hooft parameter, held fixed. More- 
over, forT=(l-e,l + e)x({|f|<5}n{t i >0}) for some e>0,6>0, 
(i) the expansion is uniformly valid on compact subsets of T; 

(ii) e g {x,t) extends to be complex analytic in T c — {(x,t) €C' +1 ||a; — 1| < e,\t\ <<5}; 
(iii) the expansion may be differentiated term by term in (x,t) with uniform error 
estimates as in (p} 
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The meaning of m is that for each g there is a constant, K g , depending only on T 
and g such that 



\ogTl gs (t)-n 2 e (x,t) ——^(x,?) 



- n, 2 9 



for (x,£) in a compact subset of T ■ The estimates referred to in (liiij) have a similar 
form with t% and ej(x,t) replaced by their mixed derivatives (the same derivatives 
in each term) and with a possibly different set of constants K g . 
Remark 2.2. Recently these results were extended to the case where j is odd \l r fj . 



In this case one should replace the normalized partiiton function ( 2.4 ) by its Szego 



representation in terms of eignevalues (3.7) 



To explain the topological significance of the e g (x,t) as generating functions, we 
begin with a precise definition of the objects they enumerate. A map E on a compact, 
oriented and connected surface X is a pair £ = (ivf(E), [i]) where 

• K(H) is a connected 1-complex; 

• [i] is an isotopical class of inclusions i : AT(E) — > X; 

• the complement of K(T,) in X is a disjoint union of open cells (faces); 

• the complement of the vertices in K(Y.) is a disjoint union of open segments 
(edges). 

When the genus of X is g one refers to the map as a g — map. What [lj effectively 
showed was that the partial derivatives of e g (l,t) evaluated at t — "count" a geo- 
metric quotient of a certain class of labelled g-maps. 

As a means to reduce from enumerating these labelled g-maps to enumerating 
g-maps, it is natural to try taking a geometric quotient by a "relabelling group" more 
properly referred to as a cartographic group [2J. 

This labelling has two parts; first the vertices of the same valence, £ have an 
order labelling l,...ne and second at each vertex one of the edges is distinguished. 
Given that X is oriented, this second labelling gives a unique ordering of the edges 



around each vertex. The fact that the coefficients of the free energy expansion (2.5| 



enumerate this class of labelled g-maps is a consequence of (2.5| (i) which enables 
one to evaluate a mixed partial derivative of e g in terms of the Gaussian Unitary 
Ensemble (GUE) where correlation functions of matrix coefficients all reduce to two 
point functions. (A precise description of this correspondence may be found in |14j). 
To help fix these ideas we consider the case of a j-regular g-map (i.e., every 
vertex has the same valence, j) of size m (i.e., the map has m vertices) which is the 
main interest of this paper. The cartographic group in this case is generated by the 
symmmetric group S m which permutes the vertex labels and m factors of the cyclic 
group Cj, which rotates the distinguished edge at a given vertex in the direction of 
the holomorphic (counter-clockwise) orientation on X. The order of the cartographic 
group here is the same as that of the product of its factors which is rn\j m . On the 
other hand the generating function for g-maps in this setting is given by 

e g (t j ) = e g (x = l,t = (0,...,0,t j )) (2.6) 

= EiKf K fH (2-7) 

■m>l '•' 

where kS 9 (m) = the number of labelled j-regular g-maps on m vertices. The factor 
— P^j perfectly cancels the order of the cartographic group, making this series appear 
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to indeed be the ordinary generating function for pure g-maps. However, for some 
g-maps the cartographic action may have non-trivial isotropy and this can create an 
"over-cancellation" of the labelling. This happens when a particular relabelling of a 
given map can be transformed back to the original labelling by a diffeomorphism of the 
underlying Riemann surface X. In this event the two labellings are indistinguishable 
and the diffeomorphism induces an automorphism of the underlying map. In addition, 
the element of the cartographic group giving rise to this situation is an element of the 
isotropy group of the given map. Hence, as a generating function for the geometric 



quotient, (2.6 1 is expressible as 



&>= £ TT^W(-^r (s) (2-8) 



e 9.. 

g-mapss 

E g (x,tj) = e g (x,t = (0,...,0,tj)) (2-9) 

1 



E 



g-mapss 

3g [X Oj 



( _, \m(E) /(E) 

|Aut(E)| l h) 



x 2 -^e g (xil 2 -H. 3 ) 



where 771(E) = the number of vertices of E, /(E) = the number of faces of E and 
Aut(E) = the automorphism group of the map E. We have included the x-dependent 



form, (2.9 1, of e g since that will play an important role later on and also to observe 
that this is in fact a bivariate generating function for enumerating (/-maps with a fised 
number of vertices and faces. Moreover, in this j-regular setting, one sees that the 
bivariate function is self-similar. This is a direct consequence of Euler's relation: 

2 — 2g = # vertices — # edges + # faces 

= m(E)-|m(E) + /(E). 

ffWx'W =x 2 - 2s (t j x j / 2 - 1 ) m ^\ (2.10) 

The presence of geometric factors such as a j is not uncommon in enumerative 

graph theory, a classical example being that of Erdos-Renyi graphs 23 1 . In the quan- 
tum gravity setting these factors also have a natural interpretation in terms of the 
discretization of the reduction to conformal structures via a quotient of metrics by 
the action of the diffeomorphism group. We refer to [TOl H] f° r further details on this 
attractive set of ideas. 

In [4], eo,ei and ei were explicitly computed for the case of valence J =4. We 
quote, from the same paper, the following conjecture (some notation has been changed 
to be consistent with ours): 

"It would of course be very interesting to obtain e g (ti) in closed form for any 
value of g. The method of this paper enabled us to do so up to g = 2, but works in 
the general case, although it requires an incresing amount of work. We conjecture a 
general expression of the form 

_ (l-zp) 2 ^ 1 (g) 

with P^ 9 ' a polynomial in zq, the degree of which could be obtained by a careful analysis 
of the above procedure. " 
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Here ZQ — zoit^) is equal, up to a scaling, to the generating function for the Catalan 
numbers; below it will signify Zofav) which is similarly related to the generating 



function for the higher Catalan numbers (4.24) 



Over the years there have been a number of attempts to systematically ad- 
dress this question by studying the resolvent of the random matrix and associated 
Schwinger-Dyson equations [TJ QI5] . Our methods take a different approach. 

The main purpose of this paper is to show how this conjecture can be verified and 
significantly extended. In particular, we will show that for the case of even valence, 
3=2v, 
Theorem 2.3. For g>2, 



,(z ) = CW + 



c^\v) 



4U(») 



(„-(*,- l)z )2 ff -2 

(Zp — l) r Q5g-5- r (zo) 

2g-l 

v-l 



{v-(v-l)zvf9-5 



(2.11) 
(2.12) 



r = max< 1, 



for all v>2. The top coefficient and the constant term are respectively given by 



C 3g-3 



M 



(5g- 5) (5g-3)v* ^ 



41iH^o, 



C (s) = -2(2#-3)! 



9-1 



(1 - S 2 ,g) V^ ( 2 ~ 2fc) 2 g_2fc + 2 



(2.g + 2)! (2g)!l2 (2g-l)l ^ {2g-2k + 2)\ 



%g) V" 

1)! ^ 



c (k) 



(2.13) 



(2.14) 



where a^ ' _ 1 (t/) (see Thm 



4- 



is proportional to the g coefficient in the asymptotic 
equation in the Painleve I hierarchy \lty and (r) m — 



expansion at infinity of the v 

r(r— 1) ...(r — m + 1). 

Our methods can be extended to the case of j odd and the derivation of the analogue 



to Theorem 2.3 is in progress (see section R3j). 

The route to getting these results passes through nonlinear PDE, in particular a 
class of nonlinear evolution equations known as conservation laws which come from 
studying scaling limits of the recursion operators for orthogonal polynomials whose 
weights match those of the matrix models. 

This appeal to orthogonal polynomials also motivated the approaches of [4] and 
[llj . However to give a rigorous and effective treatment to the problem of finding 
closed form expressions for the coefficients of the asymptotic free energy, (2.5), re- 



quires essential use of Riemann-Hilbert analysis on the Riemann Hilbert Problem for 
orthogonal polynomials that was introduced in [20j . Though we will not review this 
analysis here, we will state the consequences of it needed for our applications and 
reference their sources. 

In section [3] we present the necessary background on orthogonal polynomials and 
introduce the main equations governing their recurrences operators: the difference 
string equations and the Toda lattice equations. In section R] we describe how (2.5) 
can be used to derive and solve (in the case of even valence) the continuum limits of 
these equations which relates to the nonlinear evolution equations alluded to earlier. 
In section [5] we outline the proof Theorem |2.3| and in section [6] we describe the ex- 
tension of this program to the case of odd valence and briefly mention what has been 
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accomplished in that case thus far. This will also help to illuminate the full picture 
behind the idea of conservation laws for random matrices. 

3. The Role of Orthogonal Polynomials and their Asymptotics 

Let us recall the classical relation between orthogonal polynomials and the space of 
square-integrablc functions on the real line, R, with respect to exponentially weighted 
measures. In particular, we want to focus attention on weights that correspond to 
the random matrix weights, V(A), ( |2.2[ ) with j even. (Recently this relation has been 
extended to the cases of j odd Il7j . with the orthogonal polynomials generalised 
to the class of so-called non-Hermitean orthogonal polynomials; however, for this 
exposition we will stick primarily with the even case.) To that end we consider the 

Hilbert space iJ = L 2 (K,e _9s V W) of weighted square integrable functions. This 

space has a natural polynomial basis, {7r„(A)}, determined by the conditions that 

7r„(A) = X n + lower order terms 

7r n (A)7r m (A)e-^ V(A) dA = 0forn^TO. 

For the construction of this basis and related details we refer the reader to j8] . 

With respect to this basis, the operator of multiplication by A is representable as 
a semi-infinite tri-diagonal matrix, 



C-- 



fa 1 
b\ a x 1 


\ 


b\ a 2 





(3.1) 



V 



/ 



C is commonly referred as the recursion operator for the orthogonal polynomials and 
its entries as recursion coefficients. (When V is an even potential, it follows from sym- 
metry that Oj = for all j.) We remark that often a basis of orthonormal, rather than 
monic orthogonal, polynomials is used to make this representation. In that case the 



analogue of (3.1) is a symmetric tri-diagonal matrix. As long as the coefficients {&„} 



do not vanish, these two matrix representations can be related through conjugation 
by a semi-infinite diagonal matrix of the form diag( 1,6^ ,(6162) ,(fri&2&3) ;■■■)■ 

Similarly, the operator of differentiation with respect to A, which is densely defined 
on H , has a semi-infinite matrix representation, 2?, that can be expressed in terms of 
£ as 



v= — (c+w- 1 ) 

9s 



(3.2) 



where the "minus" subscript denotes projection onto the strictly lower part of the 
matrix. 

From the canonical (Heisenberg) relation on H, one sees that 

[9 A ,A] = 1, 

where here A in the bracket and 1 on the right hand side are regarded as multipli- 
cation operators. With respect to the basis of orthogonal polynomials this may be 
re-expressed as 



£,(£+«£•*- 



9sl- 



(3.3) 
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The relations implicit in (3.3| have been referred to as string equations in the physics 
literature. In fact the relations that one has, row by row, in (3.3) are actually suc- 
cessive differences of consecutive string equations in the usual sense. However, by 
continuing back to the first row one may recursively de-couple these differences to get 
the usual equations. To make this distinction clear we will refer to the row by row 



equations that one has directly from (3.3 1 as difference string equations. 



C depends smoothly on the coupling parameter tj in the weight V(A) (see |2.2[ ). 
The explicit dependence can be determined from the fact that multiplication by A 
commutes with differentiation by tj. This yields our second fundamental relation on 
the recurrence coefficients, 



d r 



£')_,£ 



(3.4) 



which is equivalent to the j th equation of the semi-infinite Toda Lattice hierarchy. 
The Toda equations for j = 1 are 



.9 



-g. 






b 2 

9s ri,g 3 



~b n ,g s \ a n,g s 



■fln-1, 



(3.5) 
(3.6) 



3.1. Hirota Equations 

One may apply standard methods of orthogonal polynomial theory |36j to deduce 
the existence of a semi-infinite lower unipotent matrix A such that 

C = A~ 1 eA 



where 



/0 1 

00 1 







V 



■7 



(For a description of the construction of such a unipotent matrix we refer to Propo- 
sition I of [13].) 

This is related to the Hankel matrix 



H = 



/mo mi m 2 ...\ 

mi rri2 m^ . . 
m,2 m.3 m^ . . 



V 



■7 



where 



m k - 



= / \ k e-^ v ^d\ 



is the k moment of the measure, by 



ADA^=U 

Z? = diag{do,di . 
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itll 



detH 



n+l 



defH„ 



where T-L n denotes the n x n principal sub-matrix of % whose determinant may be 
expressed as (see Szego's classical text [36]). 



detU n =n\Z {n) {t u t 2v ) 



Z^{t 1 ,t 2v ) = 



exp< 



7 J V(X m ;ti, t 2v )- 



jr- 



■cTX, (3.7) 



5 s 2^1og|A m -A £ 

where y(A;ii, t 2v+1 ) = \\ 2 + t l \+ t ^\ 2u . We set det% = 1. 

Remark 3.1. We sometimes need to extend the domain of the tau functions to 
include other parameters, such as t\, as we have done here. Doing this presents no 
difficulties in the prior constructions. 

The diagonal elements may in fact be expressed as 



'n+l,g a 



n,g B 



d n (0) 



where 



r 2 = ZM(t 1 ,t 2 v) 
Tn ' 9s Z(«)(0,0) 

Z(")(0,0) 



(3.8) 
(3.9) 



which agrees with the definition of the tau function given in ( |2.4[). T he second equality 
follows by reducing the unitarily invariant matrix integrals in (3.9) to their diagonal- 
izations which yields (3.8) [T%|. Tracing through these connections, from C to D, one 
may derive the fundamental identity relating the random matrix partition function 
to the recurrence coefficients, 



b 2 = 
n > 9 ° d 



d„ 



2 2 

T n+l,g s T n-l,g s ,2 



,(°) 



(3.10) 



which is the basis for our analysis of continuum limits in the next section. (Note that 

6q (0)=0 and therefore 6q =0.) We will also need a differential version of this 

relation: 

Lemma 3.2. (Hirota) 



<,. 



"3s 77- log 
Ot\ 

d 2 



2 
r n+l,g, 



'n,g 3 



9 , 

oil 



Z^+ l \t u t 2u ) 
Z^){h,t 2v ) 



d 2 



-'" \ogTl g3 =gl w \ogZ^\t x ,t 2v ), 



~~ 9s dt 2 



(3.11) 
(3.12) 
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(A derivation of this lemma may be found in [3].) It follows from (3.12| and (2.5| 

that 

Corollary 3.3. 



h z n%gB {x;t^ v )=X z ( z (x;t 2l ,)-\ \- ^ g z g (x;t 2 ^) + 



dt 2 



z g (x\t 2v ) — ^2 e fl (x;ii,t2i/)|ti=o 



(3.13) 
(3.14) 



is a uniformly valid asymptotic asymptotic expansion in the sense of (2.5 Hi). 

3.1.1. Path Weights and Recurrence Coefficients In order to effectively 



utilize the relations (3.3 3.4 1 it will be essential to keep track of how the matrix 



entries of powers of the recurrence operator, C? , depend on the original recurrence 
coefficients. That is best done via the combinatorics of weighted walks on the index 
lattice of the orthogonal polynomials. For the case of even potentials, the relevant 
walks are Dyck paths which are walks, P, on Z which, at each step, can either increase 
by 1 or decrease by 1. Set 



? 3 (mi,m2) = the set of all Dyck paths of length j from mi to iri2- 



(3.15) 



Then step weights, path weights and the (m 1 ,m 2 )-entry of C° are, respectively, given 

by 



w(p) = 



1 if the p step moves from n to n + 1 on the lattice 



if the p th step moves from n to n — 1 

W (p)= n w w 

steps pep 

E ^ p y 



*"mi ,7712 



(3.16) 



3.2. Dyck Representation of the Difference String equations 



The difference string equations are given (for the 2i^-valent case) by (3.3) 

£,(£ + t£ 2 »- 1 )_]=g s I. 



(3.17) 



By parity considerations, when the potential V is even, the only non-tautological 
equations come from the diagonal entries of (|3. 18 ): the (n,n) entry gives 



{c+tc 2u - r 

In terms of Dyck paths this becomes 



q s = [L + tL" *) , -(C + tC 2 "- 1 ) 

/n+l,7i \ / 71,71 — 1 



(3.18) 



e n* 



^-^+1 ^77+^ £_^ I 11 "n+4n(P)+l 11 "n+t ,„[/■'] 

P6'P 2 *'- 1 (1,0) \m=l 777=1 



i-n« 



(3.19) 



where £ m (P) denotes the lattice location of the path P after the m th downstep and 
we have used the relation x = ng s on the left hand side of the equation. 
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We illustrate this more concretely for the case of j = 2^ = 4. Referring to (3.15) 
the relevant path classes here are 

7 ,1 (1,0) = a descent by one step 

■p 3 (l,0) =paths with exactly one upstep and two downsteps (Fig. 3.1 ) 




Fig. 3.1. Elements ofV 3 (1,0) 



Note that the structure of the path classes does not actually depend upon n. This 
is a reflection of the underlying spatial homogeneity of these equations. Thus, for the 
purpose of describing the path classes, one can translate n to 0. 



Now applying (3.16) the difference string equation becomes, for n>0, 
I = (b 2 n+ i - b 2 n ) + 1 { bl +l (b 2 n + bl +1 + b 2 n+2 ) - bl (bl_ x + bl + bl +1 ) } 



where, for this example, we have set the parameter x equal to 1. 

3.3. Dyck Representation of the Toda equations 

We now pass to a more explicit form of of the Toda equations (3.4) in the case 
j = 2v: 



n dt 2v 



(£ 2 ») „ ,-{C 2v ) r 

\ / n+l,n— 1 v / Tb.n— U 



r-U b n 



7 / [ H b n+e m (P) + l 11 V;',„ i Pi 

P£V 2 »(2,0) \m=l m=l 



Once agian we illustrate these equations in the tetravalent case {y = 2). The relevant 
path class is: 

"P 4 (2,0) = paths with exactly one upstep and three downsteps. 



Applying (3.16), the tetravalent Toda equations become 



Idbl 



-^^^^(tl+in+^l+Q-intl^ + tl+^ + ^l) (3-20) 



dt 



where we have again used the relation x = ng s and then set the parameter x = l. 
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4. Continuum Limits The continuum limits of the difference string and Toda 
equations will be described in terms of certain scalings of the independent variables, 
both discrete and continuous. As indicated at the outset, the positive parameter g s 
sets the scale for the potential in the random matrix partition function and is taken 
to be small. The discrete variable n labels the lattice position on Z-° that marks, for 
instance, the n th orthogonal polynomial and recurrence coefficients. We also always 
take n to be large and in fact to be of the same order as — ; i.e., as n and g s tend to 
oo and respectively, they do so in such a way that their product 



x = g s n 



(4.1) 



remains fixed at a value close to 1. 

In addition to the global or absolute lattice variable n, we also introduce a local or 
relative lattice variable denoted by A:. It varies over integers but will always be taken 
to be small in comparison to n and independent of n. The Dyck lattice paths naturally 
introduce the composite discrete variable n + k into the formulation of the difference 
string and Toda equations which we think of as a small discrete variation around a 
large value of n. The spatial homogeneity of those equations manifests itself in their 
all having the same form, independent of what n is, while k in those equations varies 
over {—v— 1,...,— 1,0,1,..., v + 1}, the bandwidth of the (2v) th Toda/Difference String 
equations. Taking v + 1 << n will insure the necessary separation of scales between k 
and n. We define 



w = (n + k)g s 
= x + g s k = x[ 1 



(4.2) 
(4.3) 



as a spatial variation close to x which will serve as a continuous analogue of the lattice 
location along a Dyck path relative to the starting location of the path. 
We also introduce the self-similar scalings: 



S\=X 2 t\ 



S 2 „ = X 



2v 



(4.4) 
(4.5) 



that are natural given (2.101. In terms of these scalings, (3.13) may be rewritten [IB] 
as 



n,g. 



(S2u) = X [ Zq(S2 v ) H h —^Z g (s2u) + 



(4.6) 



Zg{s2u) = -T^e s (si,S2i/)| Sl =0 



and, by [H] App. A, 



(4.7) 
(4.8) 



hi 



n-\-k,g s 



(s 2v ) = x ( f (s2v,w) H h -^ g f g (s 2 ^,w) + 



fg{s2v,w)=W l 2g Z g (s 2v W U l ). 



(4.9) 
(4.10) 
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Remark 4.1. We mention here that the variables Sj as defined above differ slightly 
from their usage in related works I16[\10j where Sj =—a.jtj for appropriate parameters 
aj >0. 

We also introduce a shorthand notation to denote the expansion of the coefficients 
of /(si,S2i/+i,if) around w = x. 
Definition 4.2. For w = x + g s k, 

/( S1 , S2 ^)=£^^(^y (4.1D 

j=o ■>' v J 

where the subscript w^> denotes the operation of taking the j th derivative with respect 
to w of each coefficient of f: 



i 

g>0 



^^T.^J^u^w)^ 



As valid asymptotic expansions these representations denote the asymptotic series 
whose successive terms are gotten by collecting all terms with a common power of 1/n 



in (4-11) 



In what follows we will frequently abuse notation and drop the evaluation at 
w = x. In particular, we will write 

b u 9 ,-f:^(^y=±^£jf^^^ k y ( 4 - 12 ) 

3=0 J ' j=0 J 'g>0 

In doing this these series must now be regarded as formal but whose orders are still 
defined by collecting all terms in 1/n and g s of a common order. (Recall that g s ~ — 
so that n~ a g@ =0{n~( a+ ^)). They will be substituted into the difference string and 
the Toda equations to derive the respective continuum equations. At any point in this 
process, if one evaluates these expressions dX w — x and g s — ^ one may recover valid 
asymptotic expansions in which the a n +k,g 3 and 6„ , k have their original significance 
as valid asymptotic expansions of the recursion coefficients. 

4.1. The Continuum Limit of the Toda equations We are now in a 



position to study the Toda lattice equations ( 3.20 ) expanded on the formal asymptotic 



series (4.12) 



P£P 2 "(1,-1) \in=lj=0 J ' m=lj=0 J ' 

From now on we will take si=0, since its role in determining the structure of the 
asymptotic expansions of the b n+ k is now completed, and set s-i v = s. 

Collecting terms in these equations order by order in orders of 1/n we will have 
a hierarchy of equations that, in principle, allows one to recursively determine the 



coefficients of (4.6). We will refer to this hierarchy as the Continuum Toda equations. 
(Note that one has such a hierarchy for each value of v.) Of course this is a standard 
procedure in perturbation theory. The equations we will derive are pdes in the form 
of evolution equations in which w, now regarded as a continuous variable, is the 
independent spatial variable and si v is the temporal variable. One must still determine, 
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at each level of the hierarchy, which solution of the ode is the one that corresponds to 



the expressions given for f g in (4.10 ). This amounts to a kind of solvability condition. 
This process was carried out fully in [T^ and [12]. We will now state the results of 
that analysis. 

Theorem 4.3. 1161/ The continuum limit, to all orders, of the Toda Lattice equa- 
tions as n — > oo is given by the following infinite order partial differential equation for 
f(s,w): 

~J~ = \9s] JlfwT ■ ■ J Jwttl )• • • ) 

= / ,9s B Fg (/Ji»i/u)P)i'" ;/nj(2g+i)) 
g>0 

— C-VJ Jw i 9 S ^ 1 \J i JV) i Jww j Jwww ) ~r* " V^ - -^/ 

for(s,w) near (0,1) and initial data given byf(0,w) =w. (4-14) 

j(".s) / f \ rj(X) 

F <„> = E TTvw r " (A)+i n(¥) ( 4 - 15 ) 

where A= (Ai, A2,. . .) is a partition, with Ai > A2 > A3 > • • •, 0/ 2g + l; r,(A) = 
#{Ai|Ai=j}; ^(A) = ^ 7 r j(A) is i/ie length o/A; and |A|=^ j; Ai is the size o/A; and 



(i^' 9 are coefficients to be described in the next proposition. By (4-13), da' 

The above result effectively reduces the determination of the hierarchy to an 
enumeration in terms of a pair of partitions. The first class of partitions is fairly 
straightforward and amounts to keeping track of the partial derivatives that enter 
into the expressions at a given level. Thus at 0(g^ 9+1 ) the terms are products of 
partial derivatives of various orders (the parts) which must add up to 2g + l. These 
terms correspond to the tableaux at the (2g + l) st level of the Hasse- Young graph 
shown in the left panel of Figure |4~l] (In this we ignore, for now, the powers of 1/n 
that are internal to the asymptotic series / and its w-derivatives.) The other type 
of partition relates to how the Dyck paths enter our equations. We have already 
seen that a Dyck path is completely determined by specifying when its downsteps 
occur. The Toda equations depend explicitly on where these downsteps (the £ m (P) 
appearing in the equations as stated at the start of this subsection). However these 
two specifications can be related and the downstep times are encoded in terms of 
a partition that measures deviation of the path from a standard zig-zag path. The 
following proposition gives an explicit closed form expression for the coefficients d^' 9 
in terms of both classes of partitions. 
Proposition 4.4. J2IJ/ fApp. k.Z) 

4 = X! 2mA(/ti-»?i,...,^ + i-i?v+i) 

(i/ + 1,i/,...,2,1)C^C(2i/,2z/-1,...,i/) 

where 1Z is the set of restricted partitions (meaning that /ii > fi2 > • ■ • >H>v+l)> 
(771,. .. ,r)v+i) = {2v,2v— 2,. .. ,2,0), and m\(x\, ...,x u +i) is the monomial symmetric 
polynomial associated to A [29]. The relation of inclusion between partitions, pQfi 
means that fj,j > pj for all j . 

The right panel of Figure |4.1| exemplifies, in the case when v = A, a geometric 
realization of the partitions /1 being summed over in the above formula. Such a 



14 



CONSERVATION LAWS OF RANDOM MATRIX THEORY 



partition corresponds to a zig-zag path contained in the v x v + 1 rectangle which 
starts at the leftmost corner with the th step being a up-step and terminates at the 
rightmost corner with the (2v) th step. The red path corresponds to the distinguished 
partition r\ which records the downsteps at times 2v, 2v — 2,..., so that in this case 
the initial (0 th ) step is a downstep. The green path illustrates a typical /i which, in 
this case, takes downsteps at times (7,6,5,3,2). Given such a partition-path, one may 
project the initial point of each of its downsteps to the horizontal axis as indicated 
in the figure for the red and green paths. Then fii — rji equals the signed separation 
between the i th green point and the i th red point, reading from right to left. The 
symmetric polynomial m\ then gets evaluated at these separation values in the above 
This description could also have been formulated in terms of 



A 



formula for d 

the border strips associated to the skew tableaux /z/(i/ + l,i/, 

elaborate on that here. 



,2,1) but we will not 



\ / 



WW M' 

ff? ffl 



\ / 



m \/ H 

□ 

I 






• • • • • 

-• • • •- 



Fig. 4.1. left: A Hasse-Young graph (courtesy D. Eppstein); right: fi partition walks 



One is now in a position to deduce the form of the Toda hierarchy. This is done 
by setting i=l so that g s = l/n. One then collects all terms of order n~ 2g in the 



resulting expansion of (4.131 and this will be a partial differential equation in s and 
w that we refer to as the g th equation in the continuum Toda hierarchy. 

At leading order in the hierarchy one observes that, for general v, the continuum 
Toda equation is an inviscid Burgers equation [5S] 



~rfo — — —rd w (f a y 

ds v + l 



(4.16) 



with initial data fo—w. A solution exists and is unique for sufficiently small values of 
s. It may be explicitly calculated by the method of characteristics, also known as the 
hodograph method in the version we now present. Consider the (hodograph) relation 
among the independent variables (s,w,fo), 



w = fo + c v sfo. 



(4.17) 



Lemma 4.5. A local solution of (4-16) is implicitly defined by (4-. 11) 
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Proof. The annihilator of the differential of (4.17) 

(l + vcvsfo' 1 ) dfo -dw + c v f^ds 

is a two-dimensional distribution locally on the space (s,w,fo). An initial curve over 
the w-axis (parametrized as the graph of a function fo(w)), transverse to the locus 
where 1 + i^c^s/q - =0 locally determines a unique integral surface foliated by the 
integral curves of the vector field of the characteristic vector field 

f°=0 (4.18) 

ds 

<]W =c v f» (4.19) 



ds 

/o(0,w) = /o(w). 



(4.20) 



Equation (4.18) requires that along an integral curve of the characteristic vector field, 

dfo. 



/o is constant; i.e., 



o = ^(sMs)). 



dfo dfp dw 
ds dw ds 



dfo 

ds 



"<Wo 



udfo 



dw 



by (4.19) which is equivalent to (4.16). Using (4.20) to set fo(0,w) —w pins down our 
solution uniquely. □ 

Remark 4.6. We note that the numerical coefficients appearing in these Burgers 
equations depend only on the total number of Dyck paths in 7-' 2l/ (2,0). 
One finds [IB] from ( |4.17 ) and the self similar form of /o , 

fo(s,w) — wz (sw u ^ 1 ), that 
z o ( s ) = E ^Cj sJ where 

j>0 

'2i/-f\ . J 2v 
= (i/ + l) 



(4.21) 
(4.22) 



c„ = 2v 

c- 1 



J \3 



V- 



1 



•1 



i' + l 



and 



(4.23) 

(4.24) 



When v — 2, Q is the j Catalan number. For general v these are the higher Catalan 
numbers which play a role in a wide variety of enumerative combinatorial problems. 

4.2. Continuum Limits of the Difference String Equations The Contin- 
uum Difference String hierarchies may be derived from the difference string equations 



(3.19) in a manner completely analogous to what was done with the Toda equations 



in the previous subsection. 



Expanding ( 3.19 ) on the asymptotic series ( 4.12 ) we arrive at the following asymp- 
totic equations. 



1= E 



-2vs 



fwU) /1_ 

j\ \n 



E IIE 

Pe-P 2 "- 1 ^,-!) \m=l 3 =0 



f w w f£ m (P) + i 



iit^MP)' 



m—l j—0 
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The equations at leading order, 0{n~ l ), are 

l=d w f + 2vs ^2 vfo~ 1 dwfo 
Pev 2v - 1 {i,o) 
'2z/-l N 



or, equivalently, 



= d w f + 2i>[ ^ )sffo d w f 



d w (w - /o - c u sfo ) = 



(4.25) 



which one directly recognizes as the spatial derivative of the hodograph solution (4.17). 
Evaluating that solution at w = 1 yields 



C u SZq+Zq-1 = 



(4.26) 



which is the functional equation for the generating function of the v th higher Catalan 
numbers, mentioned in the previous subsection. 

The terms of the equations at C(n _29_1 ) can be computed directly and are found 
to have the form 



d w [f g + c v suf^- 1 f g ]+2s 



C u d w 2_^ fki'"fk v 

0<kj< g 



(4.27) 



/ 



fl-i 



fc=0 



\ fen \-k v =g 

^ „, i ^ ! +2^(^- 1) [2 5 -2]+^- 1) [2.g-4] + ...+F(- 1 )[0])=0, 



_Q w Y^ fl">>(.*i-*» j_o, lo f jp(u-l) ro „ _ 01 j_ „(u-l) 



(4.28) 



where Fg [2m] denotes the coefficient of n 2m in 

F {v-i) = V -^ r -i(A)TT/ XuKjj \ 

s Z^ FT r fAV -l-^ \ 7 1 / ' 

A:|A|=2g+13<'(A)<i/ 11 J J V >' J V 7 

4" 1,9) = X! 2m x {^i-f]i,...,^-fi v ) 

(i/,...,2 > 1)C/*C(2i/-1,...,i/) 

(% J ...,^ I/ ) = (2i/-l ) 2i/-3,...,l) 

In fact this equation is also w-exact! To see this we introduce indeterminates Xi 
and define 



f,l)( s >w) = fwU)(s,Xiw). 



(4.29) 



We also recall that the monomial symmetric function associated to the partition A is 
given by 



m\(xi,...,x l/ ) = x 1 1 ■■■x v " + 



(4.30) 
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where + . . . on the right hand side denotes adding the (minimal) number of terms 
which are gotten by permuting the variables Xi in the first term such that the resulting 
polynomial is symmetric with respect to S v . Finally, define 



f("-i) = z 
9 2g + l 



(4.31) 



E 
C(2i 



FT r(XV 
(i/,... > 2,l)C / iC(2i/-l > ... J i/-l) |A|=2. 9 llj ^ ^ eS ' 



E 

\|=2 
*(A)<i> 



En 



f (<r(i)) 
I A CT (;)! 



The following proposition shows that the O (n 2k 2ff ) terms of Fjf ' are an anti- 
derivative for the terms F¥~ [2.9 — 2k] appearing in (4.27). 
Proposition 4.7. 



Proof. By symmetry under the action of the symmetric group, it will suffice to 
check this identity for the single monomial term corresponding to x x x ■■■ ,#£". We 
directly calculate 



A/ 



« 



f (0 



(*') 



A,! ±^V ' \*f, A,! / (A^+l)! 



^n^r 1 =E(^+i)(n , 



J (X) \ ■' (X.,+1 



x X t Jw'M) j x h>+iJwV*+V_ 



i' = l 



<i^i' 



A,! 



(A*/ + l)l 



where in the second line we have made changes of variables replacing XiW by w. Using 



this we observe that the inner summands of d w Fg have the form 



a v fW i v I 

v — --^n^= v — L ^E(^'+ 1 ) IT 



|A|=2. 9 
*(A)<i> 



|A|=2.9 
£(A)<i/ 

e ^ ^ 1 w E^(^^••^^^ / ^ , 

|A| =2.9 + 1 J ' ^ l=1 

£(A)<i/ + l 






A,-! 



: (2.9 + 1) E 



1 



£(A)<^ + 1 



Ai A„ 

"*^1 ' ' ' 1*^1/ 



n 

3>l 



/«>«) 



r,-(A) 



In the second line the two inner summations of the first line have been rewritten in 
terms of clusters of partitions adjacent to a given partition of size 2<? + l. To get the 
last line we use the fact that X^> i 3 r jW = |A| =2.9 + 1. The proposition now follows 
from this observation, (4.31) and (4.28). D 
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As a consequence of this result one sees that the continuum difference string 
equation is directly integrable: 



f B = 



-2s 



l+c v svfo~ 



{ I \ 

u 2^ A-i • • • Ik 

0<kj<g 
\ fcH Vk v =g ) 



+ ,(i?r i) [2ff-2]+i?r i )[2. g -4] + ... + F(- 1 )[0])+l£ 1 - 



lkw( 2 S-2h) 



2s^ Q (2g-2k + iy. 



Setting w = l and applying (4.261 to eliminate s this reduces to 

u \ 



2z {z -l) 
(y-(v-l)za) 



E 

< kj < g 

k \ki-\ Yk v =g 



gfa 



Zk v 

zo 






/ 



Cy^Q 



?("-!) 



- Fr i; [2 5 -2]+Fr i; [2«?-4] 



("-i)r 



-^- 1} [o] 



It is immediate from this representation that z g is a rational function of zq. Apriori 
this anti-derivative should also include a constant term (in w; it could depend on 
s). This would lead to a term of the form c{s)/ (y— (v— 1)zq). However, in p~2] it is 
shown, by an independent argument, that the pole order in z a t v l^ v ~ 1) is always 
greater than one. Hence the constant of integration must be zero. With further effort 
this can be refined to 
Theorem 4.8. m 



Zg(zo)- 



zoOo-1) -£39-2(20) 



: Zq 




(4.32) 



a!f\y) 



og-iM 



(^-(i/-l)z ) 2 9 (l/-(l/-l)2 ) 2s+1 



(l/-(y- 1)^)60-1 r 



where Pz g -2 is a polynomial of degree 3g — 2 in zq whose coefficients are rational 
functions of v over the rational numbers Q and o,^ 1 - 1 (iy)^0. 

Remark 4.9. A key element in the proof of Proposition ^ -7| is the observation that 
differentiation with respect to w adjusts the multinomial labelling of partial derivatives 
in the expansion according to the edges of the Hasse-Young graph (Fig J^.l). This 
graph describes the adjacency relations between Young diagrams of differing sizes. 
The edges describe which partitions of size 2g + l are covered by a given partition of 
size 2g. Conversely it describes which partitions of size 2g cover a partition of size 
2g + 1 which in the setting described here acts as an anti- differentiation operator. This 
kind of structure was called a differential poset by Stanley and systematically examined 
in 
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5. Determining e g 

Recalling the basic identity (3.10) 



bl = r_^L lm ^ 



(5.1) 



we have, by taking logarithms, 



logr 2 +1 - 21ogr 2 + logT 2 _ 1 = log(#) -log(6*)(0) , 



(5.2) 



where the initial value 6^(0) —x is given by the recursion relations of the Hcrmite 
polynomials. As in |16j . we can use formula (5.2) to recursively determine e g in terms 
of solutions to the continuum equations. We use the asymptotic expansion of b 2 n which 



has the form (4.6): 



9=0 



-2g 



Note that the left hand side of equation (5.2) has the form of a centered second 
difference, Ait 2 n — A_it 2 n . It follows jTHj that this expression has an expansion for 
large n involving only even derivatives of the spatial variable w. We have, at order 



dw 2 



+ the n~ 2g terms of log 



Eg^l(s,w) 



(5.3) 




where Ef l (s,w) = w 2 ~ 2h ei l (w l/ ~ 1 s). In |12j it was shown that e g is rational in zq with 
poles located only at zq = -^r . However we will now prove the more refined result 
stated in Theorem 12.31 

The proof of this result is by induction on g. (The base case of g = 2 is established 



by direct calculation [16].) We assume that (2.11) holds for all k<g. We state 
here, without proof, some straightforward lemmas and propositions describing the 
derivatives of ( 2.11[ ) (details may be found in 112] where similar lemmas are proved 
for the z g ). Set 



E k {w,s)=w'-^e k {w v -'s) 



Lemma 5.1. 



whe 



Si 






3=0 

= d J ek/ds 3 '; s = sw u ~ 

=Q?Ji 1 ' k) M+{(^-i)i-(2*-3+p)}Q?- 1 ' fc) M 

= (2-2fc) p p>0 

= 1 

= j>p, j<0. 



0<j<P 
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Lemma 5.2. For l<k<g and j>0, 



/3fe-4+j (fc,j), n 



£=0 



(v-(v-l)z ) 2k + l +3- 



cp j \^ = [(j-l)u-(2k + e+(j-3))}cf j -\iy)+u(2k^+{j-3))4^ 1 (p) 
c f<V)= c fV). 



cf J) (i/)=0 £<0 7 £>3k-3+j 



Lemma 5.3. For \<k<g and j>0 

(j)_ ^0 y^ V^ \- 1 ) { r ) C ni-r 



1 ("-i^S £■ o-^-w*"- 1 



where m = £ + r. 

Proposition 5.1. For l<k<g andp>0 






(i/-(i/-i> r 



But in fact, by the following vanishing lemma 
Lemma 5.4. 



p m , .n. 

j=0 r=0 ^ ' 



for m = 0,l,...,p— 1, the minimal pole order of the expansion in Proposition |5.1| is 
>2fc — 2+ p. In particular the minimal pole orders coming from terms involving Ek 



on the right-hand side of (5.3) are all greater than 2k — 2 + (2g + 2 — 2k) = 2g. 



Proposition 5.2. The n- 2 9 terms of log(£™ =0 ^Lfm{s)) 

r 3 (A) 



l^n i >irj(A)!ll\(i/-(i/-l)«b) ! W" h '''" h (i/-( V -l)«o)W- 1 / 



6?/ (4-- 32). This result shows that the minimal pole order coming from the log terms 
in (5.3) is once again greater than 5Z 7 ->i2jrj =2g. 



The preceding lemmas and propositions provide explicit Laurent expansions (in 
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zq) for all terms on the right hand side of (5.3) with two exceptions: 
-hi\{s,w) = — — ^ei{w s 



dw 2 a' 



dw 2 a 
ei = — log(i^- (v-l)z ) 

g2g+2 g2g+2 

;E {s,w)= „ _ . 2 w e (w v a); 



9 w 2g+2 " 



5w 2 s+ 2 

1, (^-l) 2 , ^/ 3(i/ + l) 

eo = 2 logZo + M^TT) (z °- 1) ( v " 0_ ^^r 



(5.4) 



(5.5) 



With a small modification (5.4 1 may be brought in line with Proposition 5.1 
Proposition 5.3. 

( E i) w (p) ( s > !) = z o 2^ 



m=0 

(P-1)! 
12 



(i/-(i/-l)z ) 



m+l 



!/-l ^ 



m— 1 



(i/-(i/-l)^o) 



wzi/i c = ^Y2^- ^W other coefficients are then specified by the corresponding recur- 



sions stated in Lemmas \5.1\ - \5.3\ with k set to 1. 

A variant of the vanishing Lemma 5.4 also holds for {Ei) wip) (s,l): 

Lemma 5.5. 



E<rV)£(-i 

j=l r=0 



r 



(^(^^1,m) 



for m= l,...,p— 1. It follows that the minimal pole order of the expansion in Propo- 
siti on|5.3| is at least p and so the corresponding contribution to the minimal pole order 
of (|5.3| is >2g. 

Finally we observe that for p>3,-S-^Eq(s,w) is a rational function of /o and its 
w-derivatives, 
Proposition 5.4. 



Oir 



— — E (s,w )= „ I -7- +PW — 



™v/, 



(lie 



(^-l) 2 

2i/(i^ + l) 



+ (P-3)1 - 



'LIJ , X 

/ y I ■ I /omi(j) /oto(p-J) 
J=0 ^ 



2 V /o / tu(p-i) 

— J (WJq w (p) +PJq w (p-i) ) 



(5.6) 



1 



ir 



V-2 



Each line of the above proposition can be established directly by induction starting 
with the base case for p — 3. It then follows from Proposition 3.1 (iii) of [12] that the 
minimal pole order contributed by (5.5 1 is >2g. 
We are now in a position to outline the 
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Proof, (of Theorem 2.3| In JT2] (Theorem 1.3) it was shown that 

(zo-l)q d ( g )(z ) 



e g (zo) : 



(,/-(„- I)**,)"**) 



(5.7) 



where qd( g )(zo) denotes a polynomial of degree d(g) in zq. We first want to determine 
the relation between this degree and the pole order o(g). To this end we observe 



from Propositions 5.1 5.2 5.3 and 5.4 that the right hand side of (5.3), evaluated at 
1, is a rational function in zq which approaches a finite constant value as z$ - 



w 



>oo. 



From the form of the left hand side of (|5.3| evaluated at w 
asymptotic order (as zq 



1 one also sees that its 
oo ) is the same as that of e g . Hence, d(g) = o(g) — 1 and 
this shows that (2.11) is valid up to the determination of the minimal and maximal 
pole orders at Zq = v/{v—1). In the preceding lemmas and proposiitons we have 



seen that, for all terms on the right hand side of (5.3), the minimal pole order is > 2g. 



Furthermore, from these same representations together with Proposition 3. l(iii) of [T2"] 
one sees that, with the possible exception of the genus terms in (5.6), the maximal 
pole order of the terms on the right hand side of ( 5.3 ) is 5<? — 1. The apparent maximal 
pole order in (5.6) is 4g + 3 which exceeds the stated bound when g = 2,3. This 
maximal order comes from terms containing the factor f 0w (2 g +2) which are, specifically, 

1 f 0w V9+v (2i/ + l)(i/- 1) , {u-lf 

2 — 2^+1) fowi2f,+2} 2K^TT) ( /o/ "- <29+2 » ) 



u>— 1 



f Qw (2g+2) 

2/o 

f 0w (2g + 2) 



_ (2 I / + 1)( I /-1) (v-l) 2 a 

2i/(i/ + l) ° 2i/(i/+l) °. 



2/o 



| w= iO(z/-(^-l>o)- 



Hence the maximal pole order contributed by the genus terms is, in fact, 4g + 2 
which is < 5/7 — 1 for g>2 and <5<?— 1 for g>i. This establishes that, for g>2, 
the pole orders on the right hand side of (5.3) are bounded between 2g and 5g — 1. 



Moreover, for g > 3, the case by case checking of terms on the right hand side of ( 5.3 1 



that has been carried out in this subsection, shows that the maximal pole order is 



realized by the term in Proposition 5.2 corresponding to the partition A of 2g having 
minimal length (=1); i.e., the partition whose Young diagram is a single row. This 
implies that the residue of the maximal order pole is g\ ffljjg—i (y) which is non- vanishing 
by Theor em |4.8| Hence the maximal order pole is realized. Now, given that e g has 
the form (p.7| with d(g) —o(g) — 1, it follows from direct calculation that -^ I Eg{s 1 l) 
raises the minimum pole degree by 2 and the maximum pole degree by 4 with the 



coefficient at this order given by (2.13). This establishes (2.11) for g>2>. The cases of 



5 = 2,3 may be established separately by direct calculation (see, for example, section 
1.4.2 of HU). 

To establish (2.14) first note that by Euler's relation, 2 — 2g = m — vm + F for a 



g-map where m is the number of (2^-valent) vertices and F is the number of faces. 
Since F>1, one immediately sees that the number of vertices of such a map must 
satisfy the inequality 



m> 



2fl-l 

v-\ ' 



It follows that e g (s = 0) =0 for j<r = max^ 1, J'_ 1 >. (e g must vanish at least simply 
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at s = since since t^(s = 
expressed as 



= 0) = 1.) Via Cauchy's theorem these conditions may be re- 



1 
2~7ri 






ds 



l -I 



z vj - X q 5g - 6 {z) 



^ij z ~i {v-iv-^zfa-^iz-iy 



dz 



for j < r where in the second line we have rewritten e s as a rational function of z ( 5.7 1 

" V+"i which may be deduced from 



and employed the change of variables ^| 



the string equation (4.26). This yields a contour integral in z centerd at 1. Now one 

= 0for 



<S) 



= 1 



can see that these vanishing conditions are satisfied if and only if q£l_§{z 
j <r which in turn proves (2.14). 

Finally we turn to the determination of the constant C^ . By Proposition 
contributions to the constant term of e g come only from the first sum on the right hand 
side of ( 5.3 1 . The parts of this coming from g — and g 
(2s 



5.2 



5.3 



: 1 are, by Propositions 5.4 and 



respectively, —2] 



-l)! 



and 2 



(2g-l)! 1 



(2g+2)\ allu " (2g)[ 12 

to the constant term is determined by Lemma 



At higher genus, k<g, the contribution 
(2-2fe) 28 _ 2fc+2 c , (fc) _ He 



by (5.3) we have 



5.1 



to be -2- 



9-1 



(2g-2fc+2)! 



(9 2a)h 1a\rV>- o ^" 1 ) 1 i g ftg- 1 ) 1 1 y (2~ 2fc) 2g _ 2fc+2 (fc) 



(2ff + 2 



(2g)\ 12 



fc=2 



(2.g-2fc + 2)! 



from which (2.14) immediately follows. □ 



6. The Case of Odd Valence In the case when j is odd in the weight (2.3) for 



V , there is clearly a problem in applying the method of orthogonal polynomials as it 
was outlined in Section [3] Very recently, however, a generalization of the equilibrium 
measure (which governs the leading order behavior of the free energy associated to 
(2.1)) was developed and applied to this problem, [5]. It is based on generalizing 



to a class of complex valued non-Hermitean orthogonal polynomials on a contour in 
the complex plane other than the real axis. These extensions were motivated by new 
ideas in approximation theory related to complex Gaussian quadrature of integrals 
with high order stationary points [7]. 

But even when the issue of existence of appropriate orthogonal polynomials has 
been resolved, there are still a number of significant obstacles to deriving results like 
Theorem |2 ,3| that are not present when the valence j is even. For odd valence there is 
an additional string of recurrence coefficients, the diagonal coefficients a n of C, whose 
asymptotics needs to be analyzed. This in turn requires that the lattice paths used 
to define and analyze the Toda and difference string equations must be generalized 
to the class of Motzkin paths which can have segments where the lattice site remains 
fixed rather than always taking a step (either up or down) as was the case for Dyck 
paths. 

Nevertheless, all these constructions have been carried out in [T7] to derive the 
hierarchies of continuum Toda and difference string equations when the valence j is 
odd. 

The recurrence coefficients again have asymptotic expansions with continuum 
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representations given by 



g>0 

hg (si , s 2 „ +1 , w) = -w 1 ^ 9 x 

2gi+3=g + l 
9i>0,j>0 



win 



(6.1) 



«i J - 2 » l e ft f(««5)- 1 / J » ll («ni)) M / 2 * Jl(+1 



(6.2) 



The off-diagonal coefficients have corresponding representations which are much as 
they were in the even valence case, 



bl +k .N = f( s i, s 2^+i,w)^x^2fg(s 1 ,s 2 ^+i,w)n 2g 



g>0 



2 



f g (si,s 2 „+i,w)=w 2 29 ^r^e g (w 1/2 si,w" 1/2 s 2 v+i)- 

The coefficients in these expansions have a self-similar structure given by 
hg(si,s 2v +i 1 w) = w*~ 9 u g (siw~ 1/2 1 S2v+iw l '~2) 

fg{si,S2v+l 1 w)=W 1 ~ 29 Z g (siW~ 1/2 iS 2l ,+ lW U ~i). 

At leading order the continuum Toda equations are 

and the leading order continuum difference string equations are 

0\ ( A n A 12 \ d fh ' 



1 / V foA-u An I dw I /o 



where the coefficients of the matrix B in (6.7) are specified by 



2v 



fi — j. 



and those of the matrix A in (6.8) by 



v-l 



m=o 



2;/ 



Au=i+(2»/+i)«x;(o.^i „ .. i„ j^-m)^ 4 - 1 ^-"" 1 (6-ii) 



2/j,+ l > v — fj,— l,u — /j, 



A 12 = (2v + l)sJ2 



H=0 



2;/ 



2/j,,p — /j,— l,v — jjb- 



(6.3) 
(6.4) 



(6.5) 
(6.6) 



(6.7) 
(6.8) 

(6.9) 
(6.10) 



^V-M+iytfVo""" 1 - (6-i2) 
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Remark 6.1. The index [i appearing in the trinomial coefficients corresponds to the 
number of flat steps in the Motzkin paths giving rise to that term. 



It is straightforward to see that (6.7) may be rewritten in conservation law form 
as 



d_( h \ 3 / *! 
ds \fo + \hl) + dw \^ 2 + h ^ 1 



1/,-i J T M..; 1 ,T, I .7, ,T,. I - ( 6 - 13 ) 



where the coefficients in the flux vector are given by 

2^ + 1 \ h 2n f u-ii+\ 



li=0 



*i=F „ Ahffr^ 1 (6-14) 



* 2 =E 



li=0 



2V + 1 Vo^Vo"^ 1 - (6-15) 

2/j,+ l,v-fj,-l,v-ij, + lJ ° ° v ; 



Recently, jTS], we have determined that the equations (6.8) are in fact a differ- 
entiated form of the generalized hodograph solution of the conservation law (6.13). 
This hodograph solution is given by 

*i = fto + (2i/+l)«Bia = (6.16) 

$ 2 = /o + (2^+l)sBii=iy. (6.17) 



Analogous to what was done in Theorem |2.3| we expect to determine closed form 
expressions for all the coefficients in the topological expansion with odd weights. The 
first few of these, for the trivalent case, are [T7] 

,,, 1, f , , 1 (zo-l)(zg-6zo-3) 
eo(t 3 ) = - 2 lo g (z ) + - ^-— } , 

^3) = -^logg-f), (6.18) 

1 (z 2 -l) 3 (4z 4 -93z 2 -261) 
21 3J 960 (^-3) 5 

where zq is implicitly related to £3 by the polynomial equation 

1 = z — 72t 3 z . 

£0(^3) is m fact the generating function for a fractional generalization of the Catalan 
numbers. Its m th coefficient counts the number of connected, non-crossing symmetric 
graphs on 2m + 1 equi-distributed vertices on the unit circle [35] . 

7. Concluding Remarks 

7.1. Spectrum Nothing has been said, in this article, about the eigenvalues 
of the random matrix M although this is at the heart of the Riemann Hilbert anal- 
ysis underlying all of our results. The essential link comes through the equilibrium 



measure [331 US], or density of states, for these eigenvalues. When tj—0 in (2.2), this 
equilibrium measure reduces to the well-known Wigner semi-circle law. As tj changes 
this measure deforms; but, for tj satisfying the bounds implicit in (2.5) (i), its sup- 



port remains a single interval, [a,/3]. The edge of the spectrum a(tj),/3(tj) evolves 
dynamically with tj. In fact we have determined [18J that a and j3 are the Riemann 
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invariants of the hyperbolic system (6.7). In the case of even valence these invariants 
collapse to ±2^/z~q where Zq is the generating function upon which all the e g are built, 
as described in Theore m |2.3| So the edge of the spectrum is indeed directly related 
to the genus expansion ( |2.5| . 

7.2. Random Surfaces It follows from Corollary |3.3| and the map-theoretic 



interpretation of e g given by (2.8) that z g (tj) is a generating function for enumerating 
j-regular g-maps with, in addition, two legs. A leg is a univalent vertex; i.e., a vertex 
with just one adjacent edge, that is connected to some other vertex of the map. 
In particular Zo(t) enumerates such maps on the Riemann sphere or what are more 
commonly referred to as two-legged planar maps. In a remarkable paper, [34] . building 
on prior work of Cori and Vauquelin [5] , Schaeffer found a constructive correspondence 
between two-legged 2^-regular planar maps and 2v valent blossom trees. A 2v valent 
blossom tree is a rooted 2v valent tree, with each external vertex taken to be a leaf 
that is colored either black or white and such that each internal (non-root) vertex is 
adjacent to exactly v— 1 black leaves. This gives another interpretation of Zofav) as 
the generating function for the enumeration of blossom trees. It seems rea sona ble to 
hope that the arith met ic data implicit in the coefficients am (v) in Theorem 4.8 (resp. 



Cm (V) in Theorem 2.3) might provide a means, such as sewing rules, for constructing 



two-legged 2^-regular g-maps (resp. pure 2i/-regular g-maps) from blossom trees. 

In another direction Bouttier, Di Francesco and Guitter [5] have studied the 
combinatorics of geodesic distance for planar maps. They define the geodesic distance 
of a two-legged graph to be the minimum number of edges crossed by a continuous path 
between the two legs and study r^{t2 V ), the generating function for enumerating all 
two-legged 2v valent planar maps whose geodesic distance is < d. They find surprising 
and elegant closed form expressions for the r^iy). The statisitcs of planar maps is a 
natural stepping off point for the study of random surfaces. There has been a lot of 
recent activity in this direction by Le Gall and his collaborators related to the work 
of Schaeffer and Bouttier et al. See for example [28] , 

7.3. Enumerative Geometry of Moduli Spaces A different (from matrix 
models) representation of 2D Quantum Gravity may be given in terms of intersection 
theory on the moduli space of stable curves (Riemann surfaces), M g _ n and from this 
alternate perspective Witten conjectured that a generating function for the intersec- 
tion numbers of tautological bundles on M. g , n should be "given by" a double-scaling 



limit of the differentiated free energy (3.13) for matrix models. (A precise description 



of this double-scaling limit may be found in [12 j.) He further conjectured that this 
intersection theoretic generating function should, with repect to appropriate choices 
of parameter vatriables, satisfy the Korteweg-deVries (KdV) equation. Subsequently, 
Kontsevich [57] was able to outline a proof of Witten's KdV conjecture based on a 
combinatorial model of intersection theory on M. gi n- This model expresses tautologi- 
cal intersections in terms of sums over tri- valent graphs on a genus g Riemann surface. 
He was then able to recast this sum in terms of a special matrix integral involving 
cubic weights on which the proof of Witten's KdV conjecture is based. For a readable 
overview on the above circle of ideas we refer to 32J . 

However, the first Witten conjecture, on relating the intersection-theoretic free 



energy to the matrix model free energy, (2.5), remains open. With the results de 



scribed in this paper it may now be possible to determine if, and in precisely what 
sense, this conjecture might be true and to see if this leads to connections between 
the KdV equation in the Witten-Kontsevich model and the conservation laws given 
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by (6.7). In addition, given the recent results on matrix models with odd dominant 
weights, [51 [T7] , it may now be possible to give a rigorous treatment of Kontsevich's 
matrix integral which, up to now, has been formal. 

More recently there have been other, perhaps more natural, approaches to the 
proof of Witten's KdV conjecture, [25l [3TJ [21] , in terms of coverings of the Riemann 
sphere and Hurwitz numbers for which the generating functions specified in Theorem 
2.3| should also have a natural interpetation. 

7.4. Analytical Deformations and Critical Parameters In [T5] it was 



observed that the equilibrium measure (see 7.1) for the weight V, with j — 2v in (2.2), 
may be re-expressed as 

fi Vt (X) = z Q fi Q (X) + (l-z Q )fi 00 (X) (7.1) 

where [Iq is the equilibrium measure for 1^=1/2A 2 (the semicircle law) and /Zoo is 
the equilibrium measure for V = \ 2v ; i.e., the general measure is a linear combina- 
tion, over Zq, of two extremal monomial equilibrium measures. For z G [0,1] (which 



corresponds to ti v S [0,oo]), this combination is convex and (7.1 1 is indeed a measure 



with a single interval of support in A G R. This may be analytically continued to a 



complex zq neighborhood of [0,1] so that (7.1) remains a positive measure along an 



appropriate connected contour (" single interval" ) in the complex A-plane. For v = 2 
this continuation may be made up to a boundary curve in the complex z -plane pass- 
ing through zq = 2 (with a corresponding image in the complex £4 plane). Extension 
of this to more general values of v is in progress [15] . The mechanism for carrying 



out this continuation is to regard (7.1) as a zo-parametrized family of holomorphic 
quadratic differentials. The candidate for the measure's support is then an appropri- 
ate bounded real trajectory of the quadratic differential. Outside the boundary curve, 
the Riemann-Hilbert analysis used in this paper may be analytically deformed and 
our results extended. The boundary may be regarded as a curve of critical parameters 
for this deformation. This curve is precisely the locus where the Riemann invariants, 



that determine the edge of the spectrum (as described in 7.1) exhibit a shock. 

This scenario is reminiscent of that for the small fi-limit of the nonlinear 
Schrodinger equation [241 126) in which the analogue of our boundary curve is the 
envelope of dispersive shocks. In that setting it is the Zakharov-Shabat inverse scat- 
tering problem that shows one how to pass through the dispersive shocks and describe 
a continuation of measure- valued solutions with so-called multi-gap support. It is our 
expectation that coupling gravity to an appropriate conformal field theory (to thus 
arrive at a bona fide string theory) [3D] will play a similar role in our setting to deter- 
mine a unique continuation through the boundary curve of critical parameters to a 
unique equilibrium measure with multi-cut support. We also hope that this will help 
bring powerful methods from the study of dispersive limits of nonlinear PDE into the 
realm of random matrix theory. 
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